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We use the Hubbard Hamiltonian H on the honeycomb lattice to represent the va- 
lence bands of carbon single-wall (JV, N) nanotubes. A detailed symmetry analysis 
shows that the model allows W = pairs which we define as two-body singlet eigen- 
states of H with vanishing on-site repulsion. By means of a non-perturbative canonical 
transformation we calculate the effective interaction between the electrons of a W = 
pair added to the interacting ground state. We show that the dressed W = pair is a 
bound state for resonable parameter values away from half filling. Exact diagonaliza- 
tion results for the (1,1) nanotube confirm the expectations. For (N,N) nanotubes of 
length I, the binding energy of the pair depends strongly on the filling and decreases 
towards a small but nonzero value as / — > oo. We observe the existence of an optimal 
doping when the number of electrons per C atom is in the range 1.2-5-1.3, and the 
binding energy is of the order of 0.1 4- 1 meV. 



I. INTRODUCTION 

After the discovery of carbon nanotubes [1] the interest in such systems has been stimulated by their anoma- 
lous normal properties [2] and by the recently reported superconductivity [3] . Indeed there has been growing 
evidence of superconducting fluctuations in single-wall carbon nanotubes placed between superconducting 
contacts [4] [5] [6] up to the transition temperature of ~ 0.5 K [3]. 

A single- wall carbon nanotube (SWNT) is a graphite sheet wrapped onto a cylinder. The carbon atoms are 
arranged on the sites of a honeycomb lattice. The two primitive Bravais lattice vectors are a± = (d/2)(±l, v3), 
where d/y/3 is the nearest neighbor carbon separation, see Fig.(l). A SWNT is characterized by a pair of 
integers (N, M) which specifies the wrapping: the cylinder has the axis running perpendicular to iVa+ + Ma_ , 
so that atoms separated by iVa+ + Afa_ are identified. Only in recent years, was it possible to study the 
electronic properties of atomically resolved SWNT's; it was found [7] that they are strongly dependent on the 
integers JV and M. 

From band-structure calculations [8] [9] one predicts that the "armchair" (N, N) tubes are metals while the 
"zig-zag" (N, —N) ones [which are the same as the (N, 0) tubes] are insulators or semiconductors. However, the 
Coulomb interaction cannot be neglected [10] [11]. Due to the quasi one-dimensional structure, the SWNT's are 
believed to exhibit non-Fermi liquid behaviour [12]. The band structure of the armchair nanotubes suggests [13] 
[14] that near half-filling the low-energy effective Hamiltonian should be a four-component Luttinger model. 
This theory was developed by the perturbative Renormalization Group [13] [15] and by the bosonization 
technique [16] [17] in the weakly doped case, revealing the presence of a superconducting instability for short- 
ranged interactions. Still in a one-dimensional scheme, phonon induced pairing mechanisms have also been 
reported [18] [19]. 

On the other hand, in this exploratory paper we consider a different scenario in which the electron density is 
so far from half-filling that any one-dimensional free-fermion model with a linear spectrum fails to describe the 
non-interacting Hamiltonian. This is motivated by the expectation [6] that an increased doping by chemical 
manipulations and/or external fields could enhance the superconducting transition temperature. In recent 
years, we proposed a pairing-mechanism [20] [21] in the two-dimensional (one-band and three-bands) repulsive 
Hubbard model for the Cuprates. Here we show that the same idea extends to the case at hand. In our approach 
bound electron-pairs in SWNT are obtained by a symmetry-driven configuration interaction mechanism in 
which the transverse direction plays a crucial role. 

In this paper, we do not (yet) include phonons even if we acknowledge that their contribution could be 
relevant. However, on one hand, we wish to explore an electronic mechanism which per se leads to bound 
pairs. On the other hand, any mechanism in low-dimensional systems, like the Cuprates and nanotubes, 
must overcome somehow the problem of the repulsion between confined charges. This is an obvious difference 
compared to the traditional superconductors where pairs have hundredes of Angstrom of space to delocalize. 
The notion that pairing can arise by a purely electronic mechanism, i.e. from purely repulsive electron-electron 
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interactions, was put forth by Kohn and Luttinger long ago [22]. They suggested that for large odd values of 
the relative angular momentum two electrons could stay enough far apart from each other to take advantage 
of the Friedel oscillations of the screened Coulomb potential. In our approach, based on a 2d Hubbard model, 
the first-order Coulomb repulsion is removed by symmetry. 




d 
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FIG. 1. Illustration of the honeycomb lattice. The sites a (empty dots o) and b (full dots •) constitute the basis 
of each unit cell. We fix the origin in an a site. The (N, M) nanotubc is obtained by identifying sites separated by 
Na+ + Ma_ . 

The plan of the paper is the following. In Section II we define the Hubbard model Hamiltonian H for the 
armchair (N, N) SWNT and we introduce some useful notations. In Section III we use a general theorem to 
obtain all the two-body singlet eigenstates of H with vanishing Hubbard repulsion (W = pairs). We exploit 
the Space Group symmetry of the system to get all the W = pairs with zero total momentum. Remarkably, 
we find that their wave-function vanishes if the particles have the same x value. In Section IV we propose 
a non-perturbative canonical transformation to deal with the effective interaction between the electrons of a 
W — pair added to the many-body ground state. Since the two extra particles cannot interact directly by 
definition of W = pair, their effective interaction comes out from virtual electron-hole excitation exchange 
with the Fermi sea and in principle can be attractive. In Section V we consider the (1, 1) nanotube, with length 
I — 2d. This system has 8 sites and is the smallest nanotube-like cluster showing the superconducting W = 
pairing. Therefore, it represents a very good probe to test the pairing mechanism described in Section IV, since 
we can compare exact diagonalization results with the analytic ones. In Section VI we apply the canonical 
transformation approach to study the (N, N) nanotubes of finite length and periodic boundary conditions. 
We obtain a Cooper-like equation for the binding energy —A of the W = pair which is numerically solved 
for 2 < N < 6 and length I up to 32d; the results are then extrapolated to study the dependence of A on the 
radius of the tube and on the Fermi energy in the limit of infinite length. Finally the conclusion are drawn in 
Section VII. 



II. HUBBARD MODEL ON THE SWNT 

The four outer-shell electrons of each carbon atom form three a sp2 bonds and a resonant n bond with the 
remaining p z electron. A simple description consists of a tight-binding Hamiltonian where only the p z orbital 
is taken into account: 

H o = tJ2Yl ( C 'U C '-V + h.c) , (1) 

<r,r') a 

where cl >a (c r , a ) is the creation (annihilation) operator of an electron of spin a on the honeycomb site r, the 
sum runs over the pairs (r, r') of nearest neighour carbon atoms and t is the hopping parameter. The on-site 
Coulomb repulsion is 

w = u^h TA h r , y , (2) 
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where n r , CT = <i,aC r ,a is the number operator referred to the site r and to the spin a. The full Hamiltonian 
reads 

H = Ho + W. (3) 

Although the the point symmetry Group of the graphite layer is D 6h , the simplified Hubbard Hamiltonian in 
Eq.(3) has C§ v symmetry since the fermionic operators are taken to be even under reflection with respect to 
the plane; for a SWNT only C 2v point symmetry remains, due to the wrapping. The C 2v Group is abelian 
and its table of characters is 



c 2v 


1 


c 2 




(Ty 


Symmetry 


Ai 


1 


1 


1 


1 


•2 . '2 
x +y 


A 2 


1 


1 


-1 


-i 


x/y -y/x 


Bi 


1 


-1 


i 


-i 


X 


B 2 


1 


-1 


-1 


1 


y 



Table I. Character table ot the C 2v symmetry group. Here 1 denotes the identity, C 2 the 180 degrees rotation, u x and 
<j y the reflections with respect to the x = and y = axes respectively. The C 2v symmetry group is Abelian and hence 
has four one-dimensional irreducible raprcscntations (irreps) denoted by Ai, A 2 , B\, B 2 . In the last column simple 
basis functions for each irrep are shown. 

Below, we shall diagonalize the kinetic term Ho and introduce some useful notation. Let us write the 
position R of the cells in terms of a pair of integers (n, m) , R = na + + ma_ . As illustrated in Figure 1, we set 
the origin at an a site, and therefore any R translation will take us to another a site; to get a b site one must 
translate by R — (0, j^g)- It is convenient to rename the creation operators on the site r by distinguishing 
the a sites from the b sites: 

c t _ / < ^ ' = R 

r K )f * -R-to,^). (4) 

where the spin index is omitted for the sake of simplicity. Let us introduce the Bloch creation operators 



R 

with 



* R , v = ± (5) 



u ± (k,o) \ _ 1 ( 1 



u (Kb) I ^/WZ \ ± 



A(k) 



WjOl ( 6 ) 



and 

A(k) = (l+2e-^ ( ^ ) cos^). (7) 
The kinetic term Ho can be written in a diagonal form as 

fi^EE^^ttt (8) 

v k, o 

where 



r= (k) = : t\ 1 1 + In ,s' ( ^) + 4 cos (^f) cos (^'] . 



are the bonding (-) and antibonding (+) bands. 
The Hubbard interaction becomes 



W = 2NL J2 E ^i^^^(ki,k 2 ,k3,k4) C ^ C ^i<!i4; 4 , T ' ( io ) 

ki ,k 2 ,k 3 ,k 4 1/1,1/2 ,1/3,1/4 
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where 



t/^^s.-Jki, k 2 , k 3 , k 4 ) = U V (ki, C)u ( " 2) * (k 2 , C)u ( " 3) (k 3 , C)« ( " 4) * (k 4 , C)*o(ki + k 2 - k 3 - k 4 ) 

( = a, b (11) 

and 5o(k) is 1 if k is a reciprocal lattice vector G and zero otherwise; in the (N,N) nanotubes G = nG+ + 
mG_ , with n and m integers, and G± = ^ (1, ±^75)- For the (N, N) nanotubes of length I = Ld and periodic 
boundary conditions along the x direction the k-vectors are quantized as 

2n 2tt 
k x = j^m x , m x = 0, 1,... ,L - 1 ; k y = ^ Nd m v> m y = 0,1, . . . ,2N - 1 . (12) 

For k = (±^j7r, —^71") = k± the bonding and antibonding bands touch each other and e (k±) = linearly. 

In the next Section we show that the Hamiltonian in Eq.(3) admits two-body singlet eigenstates with no 
double occupancy on the honeycomb sites and we shall refer to them as W = pairs. W — pairs are 
therefore eigenstates of the kinetic energy operator in Eq.(8) and of the Hubbard repulsion W of Eq.(10) with 
vanishing eigenvalue of the latter. The particles forming a W = pair have no direct interaction and are the 
main candidates to achieve bound states in purely repulsive Hubbard models [20] [21] [23]. Incidentally, we 
note that such states are involved in the antiferromagnetic ground state of Hubbard and related models at 
half filling [24] [25] [26] [27]. 



III. W = PAIRS IN THE ARMCHAIR (N, N) SWNT 

We obtained [21] [28] a powerful and elegant criterion to get all the W = pairs. We can do that in terms 
of the Optimal Group Q of the Hamiltonian, that we define as a symmetry Group which is big enough to 
justify the degeneracy of the single particle energy levels. By definition, every one-body eigenstate of H can 
be classified as belonging to one of the irreducible representations (irreps) of Q. We may say that an irrep r\ 
is represented in the one-body spectrum of H if at least one of the one-body levels belongs to r\. Let £ be the 
set of the irreps of Q which are represented in the one-body spectrum of H. Let \ip) be a two-body eigenstate 
of the kinetic energy Ho with spin S z = 0. Then, it holds the 

W=0 Theorem: 

r)££& WP iv) \ip) =0 (13) 

where P^ v ' is the projection operator on the irrep n. In other terms, any nonvanishing projection of on 
an irrep not contained in £ , is an eigenstate of Ho with no double occupancy. The singlet component of this 
state is a W = pair. Conversely, any pair belonging to an irrep represented in the one-body spectrum must 
have positive W expectation value. 

The complete characterization of the symmetry of W = pairs requires the knowledge of the Optimal 
Group Q. A partial use of the theorem is possible if one does not know Q but knows a subgroup. It is then still 
granted that any pair belonging to an irrep of the subgroup not represented in the spectrum has the W = 
property. On the other hand, accidental degeneracies occur with a subgroup of the Optimal Group, because 
by mixing degenerate pairs belonging to irreps represented in the spectrum one can find W = pairs also 
there. This is illustrated by the example reported in Section V. 

Below, we shall apply the W = Theorem in the (N, N) SWNT to obtain W = pairs of zero total- 
momentum. Let |0) denote the electron vacuum. Exploiting the invariance of the Hamiltonian under transla- 
tions and C2v -operations the determinantal state T c^ k JO) yields nothing if projected onto an irrep of the 
Space Group with non-zero momentum. Direct inspection of the Bloch functions in Eq.(6) shows that the 
k = irreps represented in the one-body spectrum are B 2 in the bonding band and Ai in the antibonding 
band; hence, W = pairs may be obtained by projecting onto A2 and Bi. Let us consider the two-electron 
singlet state of vanishing momentum: 

V Cl ,c 2 (k,Ri,R 2 ) ee ^(OIc^c^^.tcU^+cU.^.J |0> = 

= -J= [u* (k, Ci) u (-k, ft) e lk '( R i- R 2> + u* (k, Ca) u* (-k, Ci) e"*^-^)] Xo , (14) 

where \o is a singlet spin function and for i = 1, 2 £j = o, b. The projection onto the irrep r\ of Civ yields 

^ ] jC2 (k,R!,R 2 ) = I X M (d)^ (l ,(2 (Ok,R!,R 2 ) , (15) 

6ec 2 „ 



4 



where X^'(^) i s the character in rj of the operation O of Ci v . While ?/>' Sl ' = 0, after some algebra one finds 
i> [ £l (k, Rj , Ra) = sin (fe^Xj - X 2 )) x 

x -^[t t *(k,Ci)^(-k,C2)e lfc « (yi - y2) -«*(k,C2)^(-k,Ci)e- lfc «^-^)] xo, (16) 

where R;,-,- = (X iy j,Yij). We can verify by direct inspection that i^ 2 \ 2 (k, Ri, R2) vanishes for Xi = X2, 
that is the two-body singlet wavefunction vanishes if the particles lie on the same annulus of the (N, N) tube. 
As a consequence t/)^ 2 ^ (k, Ri, R2) is an eigenstate of the kinetic energy Ho [with eigenvalue 2e(k))] and of 

the on-site Hubbard repulsion W with vanishing eigenvalue of the latter, that is i^ 2 \ 2 (k, Ri , R2) is a W = 
pair. 

W = pairs of non-vanishing total-momentum may be obtained in a similar way; however, in this prelimi- 
nary work we concentrate on pairs of vanishing total momentum. 



IV. CANONICAL TRANSFORMATION APPROACH TO THE PAIRING MECHANISM 

In this Section we intend to study the effective interaction among the electrons of a W = pair added to the 
n-body interacting ground state l^oin))- Since the two extra particles cannot interact directly by definition 
of W = pair, their effective interaction comes out from virtual electron-hole excitation exchange with the 
Fermi sea and in principle can be attractive. 

Many configurations contribute to the interacting (n + 2)-body ground state |$o(w + 2)} and we need a 
complete set 5 to expand it exactly; as long as it is complete, however, we can design S as we please. We 
can take the non-interacting n-body Fermi sphere |$o(n)} as our vacuum and build the complete set in terms 
of excitations over it. In the subspace with vanishing spin 2 component, the simplest states that enter the 
configuration mixing are those obtained from \§o{n)} by creating two extra electrons over it; we denote with 
\m) these states. Similarly, along with the pair m states, we introduce the 4-body a states, obtained from 
]$o(n)} by creating 2 electrons and 1 electron-hole (e-h) pair. Then <S includes the 6-body f3 states having 
2 electrons and 2 e-h pairs, and so on. We are using Greek indices for the configurations containing the 
electron-hole pairs, which here are playing largely the same role as phonons in the Cooper theory. By means 
of the complet set S we now expand the interacting ground state 

|*o(n + 2)) = ^o m |m) +^o a |a) +^0/31/3) + .... (17) 

m a (3 

and set up the Schrodinger equation 

tf|*o(n + 2)) = S(n + 2)|* (n + 2)). (18) 

We stress that Eq.(17) is configuration interaction, not a perturbative expansion. When the number n of 
electrons in the system is such that |$o(^)} is a single non-degenerate determinant (the Fermi surface is 
totally filled), we can easily and unambiguously define and calculate the effective interaction between the two 
extra electrons since the expansion in Eq.(17) for the interacting ground state is unique: this is done by a 
canonical transformation [20], [29], [34] from the many-body Hamiltonian of Eq.(3). We consider the effects 
of the operators Ho and W on the terms of |^o(n- + 2)). Choosing the m, a, f3, ... states to be eigenstates of 
the kinetic energy Ho we have 

Ho\rn) = E m \m), H \a) = E a \a), H \/3) = E p \0), .... (19) 
Since W can create or destroy up to 2 e-h pairs, its action on an m state yields 

W\m) = Yw™',m\m') + J2 W a , m \a) + ^ Wf,, m \0). (20) 

m' a. f3 

The action of W on the a states yields 

W\a) = Y^W m , a \m) +J2 W c.\a\a) +J2 W ^) +J2 W ^' (21) 

ma' (3 7 

where scattering between 4-body states is allowed by the second term, and so on. In this way we obtain an 
algebraic system for the coefficients of the configuration interaction of Eq.(17). However to test the instability 
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of the Fermi liquid towards pairing it is sufficient to study the amplitudes a m of the in states. In the weak 
coupling limit this can be done by truncating the expansion in Eq.(17) to the a states because, as we have 
shown [29], the inclusion of the /3, 7, . . . states produces a S-dependent renormalization of the matrix elements 
of higher order in W, leaving the structure of the equations unaltered. 
By taking a linear combination of the a states in such a way that 

(H + W) ata ,=S aa ,K (22) 

the algebraic system reduces to 

[E m - E(n + 2)] a m + ^a m , VK m , m , + ^a a W m , a = (23) 

m 1 a. 



[K - E(n + 2)] a a + ^<wW a , m , = 0. (24) 

m 1 

Solving for a a and substituting in Eq.(23) we exactly decouple the 4-body states as well, ending up with an 
equation for the dressed pair \ip) = a m \m). The effective Schrodinger equation for the pair reads 

(H + W + S[E}) |V) = i? P ai r jV> = *W (25) 

where 

(S[E]) m , m ,=-J2 Wm E ?™ a E m ' - (26) 

is the scattering operator. The matrix elements W m ^ m i in Eq.(25) may be written as the sum of two terms 
rep 

W„ 



representing the direct interaction Wif} ' , among the particles forming the pair and the first-order self-energy 



W m , m , = W^ m , + S m , m ,W m . (27) 

Analogously in S[E] we may recognize two different contributions; one is the true effective interaction W e s 
between the electrons of the m states, while the other one is the forward scattering term F 

S m ,m' = (W / cff)m,m / + F m &m,m' ■ (28) 

The first-order self-energy and the forward scattering term are diagonal in the indices m and m'. W m and F m 
renormalize the non-interacting energy E m of the m states: 

E m -> E^ = E m + W m + F m . (29) 

Eq.(25) is of the form of a Schrodinger equation with eigenvalue E(n+2) for the added pair with the interaction 
W w + W cS . Here the W = pairs are special because W 7 ^' vanishes. We interpret a m as the wave function 
of the dressed pair, which is acted upon by an effective Hamiltonian Hp^. This way of looking at Eq.(25) is 
perfectly consistent, despite the presence of the many-body eigenvalue E(n + 2). Indeed, if the interaction is 
attractive and produces bound states the spectrum of Eq.(25) contains discrete states below the threshold of 
the continuum (two-electron Fermi energy) . This is a clear-cut criterion for pairing, which is exact in principle. 
The threshold is given by 

£f =im[£<?(£)], (30) 

{m} 

which contains all the pairwise interactions except those between the particles in the pair; it must be deter- 
mined once Eq.(25) has been solved (since F depends on the solution). The ground state energy E may be 
conveniently written as E^JP + A. A < indicates a Cooper-like instability of the normal Fermi liquid and 
its magnitude represents the binding energy of the pair. 

We emphasize the fact that in principle the canonical transformation is exact because in this way our 
framework does not require U /t to be small. The next problem is how to find a practical estimate of the 
renormalized Fermi energy. In the numerical calculations, some approximation is needed. In Section V and 
VI, we shall compute the bare quantities; that is, we shall neglect the 6-body and higher excitations in the 
calculation of W e fF and F. This is a resonable approximation if we compute small corrections to a Fermi liquid 
background and the exact numerical results in the (1, 1) nanotube suggest that this is the case, see Section V. 

We want to stress that once the expansion of |^o(^ + 2)) in Eq.(17) is truncated as specified above we do 
not need to construct a good approximation of the interacting ground state wave function in order to get the 
a m amplitudes at weak coupling; in this way we obtain information about pairing. 
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V. THE (1, 1) NANOTUBE: EXACT DIAGONALIZATION AND PAIRING MECHANISM 



There is evidence from cluster calculations that pairing can arise in the repulsive Hubbard model [23] [30] 
[31]. In this Section we shall consider the (1, 1) nanotube, with I = 2d and periodic boundary conditions, 
see Fig(2.a). This system has 8 sites and is the smallest nanotube-like cluster showing the superconducting 
W = pairing. Therefore, it represents a very good probe to test the pairing mechanism shown in Section 
IV, since we can compare exact diagonalization results with the analytic approximations of the canonical 
transformation. 

We define, following Refs. [32] [33], 

A(n + 2) = E(n + 2) + E(n) - 2E{n + 1). (31) 

where E(n) is the ground state energy with n electrons (referenced to the electron vacuum). |A(n + 2)| is one 
definition of the pairing energy. This definition is simple, but requires computing the eigenvalues with great 
accuracy, and has several drawbacks. It says nothing about the dynamics which leads to pairing. Moreover, 
generally a negative A does not unambiguously imply pairing, and further problems arise since the above 
definition depends on the comparison of systems with different n. 

However, in several studies of W = pairing in finite systems when it was possible to compute A by exact 
diagonalization we pointed out [21] [23] [34] that at least at weak coupling it agrees well with A as obtained 
by the canonical transformation. This supports the application of Eq.(31). 

Below we perform a group-theoretical analysis and obtain W = pairs by exploiting the W = Theorem. 
Next, we compute the interacting ground state energy with 2, 3 and 4 electrons by exact numerical diagonal- 
ization in order to get A(4). Finally the canonical transformation is applied to evaluate A(4) which will be 
compared with A (4). 

A. Symmetry Properties and W = Pairs 

The First Brillouin Zone (FBZ) consists of 4 points and since there are two atoms in the unit cell, two bands 
result. The symmetry properties of this system are intriguing: it is not only invariant under the operations 
of the Space Group (translations and (^-operations), but also under the dynamical operation d shown in 
Fig.(2.b). The dynamical operation d is reminiscent of a similar symmetry which must be taken into account 
to understand the degeneracies in the 4x4 Hubbard model [21]. 



7 5 8 4 




FIG. 2. (a) The Hubbard model of the (1,1) nanotube; the dashed lines denote hopping interactions due to the 
wrapping, (b) Illustration of the dynamical symmetry d\ the sites 1, 8, 7, 6 undergo a clockwise rotation while the sites 
3, 2, 5, 4 a counterclockwise rotation. One can see by direct inspection that the nearest neighbours of each site are the 
same of Fig. (2. a). 

Including d and closing the multiplication table we obtain a symmetry Group Q with 48 elements in 10 
classes as shown in Table II. 



Ci 


C 2 


C 3 


Ca 


C 5 


c 6 


C 7 


C 8 


Cg 


Cio 


1 


c 2 


Ox 


<7y 




t 1 

H,0 




d 


C 2 d 


a yd 



Table II. Top row: symbols of the 10 classes d of Q; bottom row: one typical operation for each of the classes; the 
others can be obtained by conjugation. The operations are: the identity 1, the translation t„. m of n steps along a+ and 
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m along a ; d is the dynamical symmetry. The other operations C 2 , C4 , c x , &y are those of the Group of the rectangle 
and are referenced to the centre c of Fig(2.a). 



In Table III we report the character table of the full symmetry Group Q of the (1, 1) nanotube of length 
I = 2d. 
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Ci 


c 2 


c 3 


c 4 


c 5 


c f> 


c- 


c 8 


c 9 


Cl() 


A! 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


A 2 


1 


1 


-1 


-1 


-1 


1 


-1 


1 


1 


-1 


Bi 


1 


-1 


1 


-1 


1 


1 


-1 


-1 


1 


1 


B 2 


1 


-1 


-1 


1 


-1 


1 


1 


-1 


1 


-1 


E 1 


2 





2 





2 


2 








-1 


-1 


E 2 


2 





-2 





-2 


2 








-1 


1 


Ti 


3 


1 


-1 


1 


3 


-1 


-1 


-1 








T 2 


3 


-1 


-1 


-1 


3 


-1 


1 


1 








T 3 


3 


-1 


1 


1 


-3 


-1 


-1 


1 








T 4 


3 


1 


1 


-1 


-3 


-1 


1 


-1 









Table III. Character table ot the Optimal Group Q. The irreps Ai, A 2 , B\ and B 2 reduce to the corresponding ones 
of the subgroup C 2v if Q is broken. 



The one-body eigenenergies for t = 1 eV are shown in Table IV together with the irreps of the associated 
eigenvectors. From Table II and III we see that Q is an Optimal Group as defined above. 

Quasi-momentum Energy (y — — ) Irrep Energy (y — +) Irrep 

ki = (0,0) -3 B 2 3 A 1 

k 2 = (0,-2§) -1 Ti 1 T 3 

k 3 = (7r,0) -1 Ti 1 T 3 

k 4 = (-7T,0) -1 3\ 1 T 3 

Table IV. One-body spectrum for i = 1 eV. The energies are in eV. 



In Section III we have shown how to get all the W = pairs of vanishing momentum and belonging to 
the irrep A 2 of C 2v . However, the (1, 1) nanotube is too small to achieve this kind of W = pair-states and 
the projection on the irrep A 2 of c£ |C^ k JO) is identically zero. Nevertheless, this cluster admits W = 0-pair 
solutions of different type and they may be obtained by applying the W = Theorem [28], which contains a 
very general prescription to determine all the W = pairs, see Section III. 

In the (1, 1) nanotube W — pairs formed by particles of the same kinetic energy may be obtained in the 
3-fold degenerate one-body levels with energies ±1 eV by projecting c^. T Cj c . JO) with i,j = 2,3,4, onto the 
irreps which are not represented in the one-body spectrum. In this way we find two singlets 

l^ Bl > = ^(cJ C 3, T 4 3 ,|-4 4 ,T4 4 ,t)|0> (32) 

\^2 El ) = ^(2 4 2 , T 4 2 ,i -4 3 ,T4 4 ,t-4 4 , T 4 3 ,t)lo> (33) 

and three triplets 

l^ T2 ) = ^(4 3 ,T4 4 , i -4 4 ,T4 3 ,i)lo) (34) 

l^ T2 >= ^(4 3 , T 4 2 a-4 2 , T 4 3 a)|0> (35) 

l^ T2 )= ^(4 4 , T 4 2 a-4 2 , T 4 4 ,i)|o> (36) 

while all the other projections yield nothing. 

In the above expressions the Bloch states c£ JO) can be taken both in the bonding and the antibonding 
bands (e = =pl eV respectively). In conclusions we found that the (1,1) nanotube with 8 sites has two W = 
pairs in each band. As already observed, such pairs have a non-vanishing total-momentum, contrarily to the 
W = pairs of Section III. 
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B. Exact diagonalization data 



We have performed the numerical diagonalization of the Hamiltonian in Eq.(3) for the (1, 1) nanotube of 
length 2d and periodic boundary conditions, filled with 2, 3 and 4 electrons. The interesting case arises when 
the number of electrons is n + 2 — 4, since the non-interacting ground state |<J>o(2)} is a non-degenerate singlet 
with vanishing momentum containing two electrons in the lowest energy level. Therefore, the two added 
electrons may form a W = pair according to the result of Section V A. 

We found that the interacting ground state |>I'o(2)} with 2 electrons is a totally symmetric singlet with 
vanishing momentum for any value of the ratio U/t. On the other hand, 1*0(3)} is three- fold degenerate in 
the sector S z = |, with momenta (0,-^=), (tt, 0) and (— 7T,0), for U/t < 10 5 . It belongs to the irrep Ti in 
the weak coupling regime. With four electrons (784 configurations) the ground state is a doubly degenerate 
singlet, with momenta (0,0) and (0, -^|), and belongs to the irrep 1 E\. This is the symmetry of the W = 
pair; our approach predicts the correct symmetry of the ground state. The first excited state with 4 electrons 
is a three-fold degenerate triplet with momenta (0,0), (n, 0) and (— 7T, 0) and belongs to T2. We show below 
that we are able to predict the symmetry of this state as well. We found that there is no level-crossings up to 
U/t < 10 5 . 

In order to study the pairing between the two added electrons, we must compute the quantity A(4) = 
£(4) + E(2) — 2E(3), which is related to the effective interaction according to the discussion made at the 
beginning of Section V. A(4) has been computed in a large range of U /t values, and its trend is shown in 
Fig.(3). 



5(4) 




Log U/t 

FIG. 3. Trend of A(4) versus Log U/t for U/t in the range 1 -MO 5 , A(4) is in cV. 

For arbitrary small values of U/t, A(4) is negative and decreases up to a characteristic value of U/t ~ 4 -7- 5, 
where a minimum is reached; at the minimum A(4) ~ — 0.018eV. Thus, as in the model Cu-O clusters [23], 
A(4) is a new energy scale of the system, widely different from any of the input parameters. For larger U/t, 
A(4) increases until A(4) = for U/t ~ 14. As far as U/t — > 00, A(4) increases monotonically up to the 
asymptotic positive value of ~ 0.063 eV. We emphasize that A(4) becomes positive for large values of U/t and 
hence pairing disappears in the strong coupling regime. Therefore, the above pairing mechanism cannot be 
related to the ones considered within the framework of the t — J- like models, where an infinite U is required to 
forbid double occupation on the same site and pairing is achieved by means of residual attractive interactions. 



C. Analytical Canonical Tranformation: Pairing Mechanism 

In this Section we study the (1, 1) nanotube with 8 sites by implementing the canonical transformation 
described in Section IV but truncated to the a-states. We shall compare the analytic results with the numerical 
data obtained previously. The non-interacting Fermi sphere |$o(2)) = ^.c^ ^ |0) has two electrons in the 
lowest level. Since the one-body levels are widely separated the intra-shell interaction is much more important 
than the inter-shell one. Therefore, we consider only the |m) states having the added pair of electrons in the 
lowest unoccupied level of energy -1 eV, neglecting the higher-energy orbitals. We recall that in the three- fold 
degenerate level with energy -1 eV we can write 5 two-body states (2 singlets and 3 triplets) without double 
occupancy, see Eqs. (32-36). Hence, we may further reduce the set of the m states in the expansion of the 
interacting ground state |*o(4)} [see Eq.(17)] by dropping the ones with non-zero direct interaction (which are 
expected to have wrong symmetries to be the ground state). A convenient basis for the analytic evaluation of 
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the matrix elements W m>m i and (S[E]) m . m i (see below) is obtained in terms of the creation operators 

/l,CT = ^/=( C ic 3 ,CT C k4,0")' /2,CT — ^=( C ic3,CT _ ^4,^)- (^) 

The determinantal states 

l"»i>=/ 1 t , T / 2 t , i |*o(2)>, \m 2 ) =fltft d |*o (2)) (38) 

have projection only on the the first component of -Bi and on the first component of T2, and are mixed by 
the operators W and S[E] in Eq.(25). Indeed, the symmetric combination of \mi) and |ra 2 ) yields \ip 1 1 ) of 
Eq.(32), while the antisymmetric combination yields \ip 1 T ' 2 ) defined in Eq.(34). Hence, the eigenvalue Eq.(25) 
reduces to: 



Ami \ / & 



E mi + W mi + F mi (W cB ) mi ,rri2 

Fm.2 J \ y I 61 



(39) 



where we have taken into account that = since \m\) and |m 2 ) have no-direct interaction and that the 

diagonal part of the effective interaction W c « vanishes due to the Pauli principle 1 . We performed the sum in 
Eq.(26) over the a states analytically using non renormalized a-state energies, which is justified at least in 
the weak-coupling regime. As a consequence of the fact that (W e s)m,m = the forward scattering term turns 
out to be given by 

F m = (SlE ])m , m = £ (± + ^J) . f[E\ (40) 

and it is independent of the pair index m. The first-order self-energy W mi is equal to W m2 and we defer the 
reader to Section VI for the proof in a more general case; here we limit to write the final result 

W F = W mi =W m2 = —. (41) 

The off-diagonal matrix elements (W e s)m 1: m 2 = (Weffjmj.mi = w e ff[-E] are responsible for a spin-flip-like 
effective interaction due to the structure of the states mi and mi. Performing the sum over the a states and 
taking the bare energies E a we get 

Hence, the eigenvectors of Eq.(39) are ^j(l, 1) and -^(1, —1) and correspond to the states 

|*l Bl > = j=(\mi) + |m 2 » = ® |* (2)>, (43) 

|*i T2 ) = j=(\mi) - |m 2 » = |^ T2 ) ® |<l>o(2)). (44) 

From the symmetric combination we get the singlet belonging to the irrep 1 Ei and the lowest energy Es 
satisfies the equation 

-8t + W F + f[E s ] +w eS [E s ] = E s . (45) 

According to the treatment described in Section IV we write Es = — St + Wf + f[Es] + A[_Es] so that 

A[E S ] = w eS [E s ] (46) 

The antisymmetric combination corresponds to the triplet belonging to the irrep T2 and the lowest energy Et 
satisfies the equation 

-8t + W F + f[E T ] - w e s[E T } =E T . (47) 

In this case — A[Et] = w e s[E T ] and it coincides with the modulus of w c g[Es] in standard second-order 
perturbation theory [23], see Fig(4.a). 



1 More generally, (W e s)m,m — if the m state is a determinantal state. Let m = (<£>i,t, ¥>2,|) be a pair index 
containing two spin-orbitals; drawing the diagram for the effective interaction W e « with incoming <pi^ and 
ifi 2 ,i and outgoing ipi^ and y> 2 ,i one realizes that one of the two spin-orbitals must be occupied twice. 
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(a) (b) 
FIG. 4. (a) -A[Es] and A[E T ] versus (7/t. We note that for U/t > 1 it is evident that -A[E S ] + A [Ex] ■ (b) 
Comparison between A[Eg] and A(4). 



In Fig. (4.6) it is reported the comparison between A(4), as denned in Eq.(31), and A[Es] obtained from 
the canonical transformation. We have exactly fitted the 8 points for A(4) and A[£?s] by using a polynomial 
of order 7 in U/t. The ratio r between the quadratic coefficient of A[Es] and the quadratic coefficient of A(4) 
is r = 1.00003, while the linear coefficients are essentially zero in both cases. We observe that the analytical 
value |A[_Es]| is ~ 2 times greater than |A(4)| for U/t ~ 1. This means that the inter-shell interactions 
and the renormalizations of the a-state energies have an important weight in determining the right value of 
A[_Es]. However, what is comfortable is that the analytical approach predicts the right trend of the binding 
energy: the singlets feel attraction, while triplet repulsion. This is the result we need to apply the canonical 
transformation to larger and more physical systems. 



VI. PAIRING IN NANOTUBE-SUPERCELLS 



A. Canonical Transformation and Pairing Mechanism for (N, N) Nanotubes 



Recently, Potassium [35] and Lithium [36] have been intercalated in single- and multi-wall carbon nanotubes 
and a net charge transfer was observed between the alkali-metals and the carbon atoms. Since the intercalation 
causes a very small structural deformation, we may say that the alkali-metal electrons fill the original bands of 
the nanotube carrying the system away from half-filling. In this Section we apply the canonical transformation 
to study the pairing mechanism in the electron-doped (sf > 0) armchair nanotubes of length / = Ld and 
periodic boundary conditions. Omitting the band index we choose the m states as determinantal eigenstates 
of the kinetic energy Ho- 



cLltcL,!!* ^))' e(ki),e(k 2 ) > e F 



(48) 



where e F is the Fermi energy and \<&a(n)) is the non interacting Fermi sphere with n particles. The first-order 
self-energy W m is given by 



w m = u 



^Kk,oi>(k',OI 2 + ^Wk,0 



2 

j^Kk^c)! 2 



= U 



n/2 



+ 



n/2 
4N 2 L 2 



W F 
(49) 



where jm) = c\, j.c£ jJ$o(>i))j the sums over k and k' run over the occupied states and we used the fact that 
|tt(k,C)| 2 = 1/(4NL) for any k and £. 

We truncate the expansion of the interacting ground state |\£ r o(rc + 2)) of Eq.(17) to the a states and a basis 
for them looks like 



t ki,T t k 2 ,.l 



f'k. 



|$o(n)), e(ki),e(k 2 ),e(k 4 ) > e F > e(k 3 ), (r =T,1 



(50) 



In this scheme of approximation, we want to obtain the effective interaction between the electrons forming 
the W — pairs studied in Section III. Projecting on the irrep A2 the vanishing momentum pair |-0(k)) = 



c k,T c -k,il $ °( n )} we g et from Eq.(26) 



<^ lA2l (k)|S[£]|V lA2l (k')) = -25(k-k') 



occ cmp 

d(e(k + p - q) - e F ) \U„[k, p, k + p - q, q) 

p,v q 



+ 



2 E* 



e(k + p + q) - e"(p) + e(q) + e(k) - E 

^(C)Ye(:(Clt' hk+p) Cj) ^(Q k/ + k + p^- k .Q k ' ; p) t/ 4k,p,ok , +k + p,-6k / ) 



e(Ok' + k + p) - e"(p) + e(k') + e(k) - E 



(51) 
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where we have taken into account Eq.(50) for the a states and we set E' a = E a , which is justified at least in 
the weak-coupling regime. The sums are over occupied p and empty q. Since ef > 0, we omitted the band 
index + of the electron eigenenergies e, while we sum over the band index v of the virtual hole. The vertex 
U v is given by Eq.(ll) where v refers to the Bloch-function depending on p, while all the other band indices 
are intended to be +. The first term on the r.h.s. of Eq.(51) is the forward scattering contribution F(k), 
while the second term represents the effective interaction W e g(\s.,\s.'). We may see that standard perturbation 
theory yields the arithmetic mean of the two unperturbed limits E — > 2e(k) and E — > 2e(k'). In particular, 
the forward scattering F in Eq.(51) coincides with the second-order self-energy of the one-particle propagator 
in the same limit. 

We emphasize that Eq.(51) is characterized by a symmetry-induced quantum mechanical interference of 
several terms. This interference produces a partial cancellation, and the absolute value of the result is typically 
much smaller than individual contributions. This means that the interaction is dinamically small for W — 
pairs: they have no direct interactions (that is W^, defined in Eq.(27), is zero for W — pairs), and because 
of the interference the effective interaction is reduced compared to what one could expect by a rough order-of- 
magnitude estimate. However, the presence of the theta functions and the anisotropy of the integrands prevent 
a total cancellation. Substituting Eqs. (49-51) into Eq.(25), we cast the result in the form of a Cooper-like 
Schrodinger equation 

[2e(k) + W F + F(k)] a k + J2 W ° s ( k ' k ') a * = Eak ■ ( 52 ) 

k'eB/4 

for a self-consistent calculation of E (since W e s and F are ^-dependent). The indices k and k' run over 1/4 
of the empty part of the FBZ and we denoted such a set of wavevectors as 25/4. We are interested in the 
possibility that E = 2ef + Wf + -F'min(kF) + A, with a positive binding energy —A of the W = pair; here 
■FminCkp) is the minimum value of F(k) among the kf- wavevectors on the Fermi surface to be determined 
self-consistently with the eigenvalue E. 

We have performed numerical estimates of A by working on supercells of 2 N x L = Nc cells. Here we 
solved the Cooper-like equation in a virtually exact way for N up to 6 and L up to 32. The results have 
been reported in Table V.a, V.b, V.c and V.d where the hopping parameter t = 1 eV. As a reasonable value, 
we have taken U/t = 1.7 which is of the correct order of magnitude for graphite [37] [38]. In line with our 
previous finding in the small (1,1) cluster, W = A2 singlets show pairing. The calculations are performed 
with the Fermi energy ef varying between 0.5 eV and 1.2 eV (half filling corresponds to ef = 0). We see that 
the binding energy — A of the pairs decreases monotonically both with the radius and the length of the tube. 

We need to work away from half filling in order to operate our mechanism; on the other hand, close to half 
filling the system is a Luttinger liquid down to extremely low temperatures where a gap could open [14] . 
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-V 


L N 
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12.0 


0.43 


10 4 


6.3 


0.37 


10 6 


4.0 


0.39 


20 


2 


5.7 


0.38 


20 4 


3.1 


0.38 


20 6 


2.4 


0.47 


32 


2 


3.6 


0.35 


32 4 


2.0 


0.35 


32 6 


1.5 


0.38 



Table V.a. Data at e F = 0.5 cV; -A is in meV; V is in eV, t = 1 eV and U = 1.7 cV. 



L N -A -V 
10 2 11.6 0.42 
20 2 6.0 0.40 
32 2 3.7 0.36 



L N —A —V 

10 4 6.0 0.40 

20 4 3.0 0.40 

32 4 2.0 0.43 



L N -A -V 

10 6 4.2 0.38 

20 6 2.2 0.36 

32 6 1.1 0.30 



Table V.6. Data at e F = 0.8 cV; -A is in meV; V is in eV, t = 1 eV and U = 1.7 cV. 



L N -A -V 
10 2 11.1 0.40 
20 2 6.0 0.40 
32 2 3.9 0.38 



L N -A -V 

10 4 6.5 0.24 

20 4 2.9 0.15 

32 4 1.8 0.14 



L N -A -V 
10 6 4.4 0.17 
20 6 3.5 0.25 
32 6 1.6 0.18 



Table V.c. Data at e F = 1.0 cV; -A is in meV; V is in eV, t = 1 eV and U = 1.7 cV. 
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20 4 
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0.34 
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3.9 


0.45 


32 4 


1.9 


0.37 


32 


6 


1.0 


0.33 



Table V.d. Data at e F = 1.2 eV; -A is in mcV; V is in cV, t = 1 eV and U = 1.7 eV. 



B. Extrapolation to large L 



With supercell sizes Nc > 400 numerical calculations become hard. Since we are concerned with the 
asymptotic behaviour for fixed N and L — > oo and — A(N, L) depends on N and L in a complicated way, we 
need a method to make reliable extrapolations of the numerical results. To this end, like in previous work [20], 
[31] we define the Average Effective Interaction V. This is such that setting in Eq.(52) W e g = — 7^7, with a 
constant V > for all k and k' in T>/4, one obtains the correct value of A. In other terms, once the binding 
energy — A is known by solving Eq.(52), the constant V must be chosen in such a way that 



V N c ^ [2e(k)+F(k)] 

ke-D/4 



[2e F +F min {k F )]-A(N,L)' 



(53) 



In Table V.a, V.b, V.c and V.d. we have reported the value of V which remains fairly stable around ~ 0.4 
eV and does not drop to in the limit of large L. Therefore V must be interpreted as a characteristic energy 
scale of the system which is largely independent on the Fermi energy and on the radius. We have numerically 
observed that |-F(k) — -F'min(k.F)! <C —A; hence we may extrapolate the asymptotic value of —A from Eq.(53) 
by dropping the difference F(k) — J ? m i n (k_F) and taking the limit L — -> oo. The domain is extended to T> and 
the result is divided by 4. Since k x is the component along the tube axis, we convert the summation into an 
integral: 



1 1 



r- > dk. 



V 8N 2tt ^ 

k 



1 (e(k x , k y ) - e F ) 



2(£"(&x, ky 



*asympt 



(AO 



(54) 



with Aasymp^iV) = UmL~*oo A(N,L). We solved the above equation for the unknown A asymp t(A r ) for several 
values of N and e F , using a typical value V ~ 0.4 eV obtained from the calculations in supercells. We found 
that A a sy m pt is strongly dependent on the filling at fixed N. Remarkably, there exists an optimal doping at 
the Fermi energy e F ~ 1 eV where A asym pt(iV) is appreciably different from zero, while A asym pt is strongly 
suppressed far from it for N > 4, see Fig. (5. a). 



A asympt 




(a) 



A asympt 



FIG. 5. Results of the canonical transformation approach with t = 1 eV and U = 1.7 eV. (a) —A 



asympt 



N 



as a function 



of the Fermi energy e F for N = 4 (empty triangles), N = 8 (black boxes) and N = 16 (grey diamonds). The Fermi 
energy varies in the range 0.6 -r 1.2. (b) — A aS ympt as a function of N for N in the range 4-^36 and ep = 1 eV. 



The existence of an optimal doping can be understood by looking at the density of states p (e) and at the 
integrand l/[2(e — e F ) — A asy mpt ] = g(e) in Eq.(54). From Fig. (6) we see that p(e) has an absolute maximum 
at £ = 1 eV, which is reminiscent of the Van Hove singularity in the graphite sheet, and exhibits well-defined 
oscillations as a function of the nanotube radius [39]. On the other hand, g(e) is peaked at e = e F . Therefore, 
for e F ~ 1 eV we have a synergy which leads to an absolute maximum of — A asympt . The prediction of 
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an optimal doping follows from Eq.(54) and is therefore largely mechanism-independent. A trend similar to 
the one shown in Fig.(5.b) was reported by Benedict et al. [40], although they considered a phonon-driven 
mechanism. 



P(e) 



P(e) 



sit 

(a) (b) 

FIG. 6. Density of states in the (4, 4) (a), 



P(e) 




AM 




(c) 



(ft) and (16, 16) (c) nanotube. p(e) is in eV . 



At the optimal doping we observe that — A aaym pt(iV) decreases monotonically as the radius \/lidN/2iY of the 
tube increases, see Fig. (5. b). However, in the limit of large N, A asym pt(JV) remains stable around 1.7 meV 
and may be interpreted as the binding energy of the W = pair in an optimally doped graphite sheet. By 
a rough order of magnitude estimate, we may say that the superconducting critical Temperature predicted 
by our approach at ef = 1 eV (which corresponds to a number of electrons per graphite atom of 1.25) is 
T c w Aasympt(oc) ~ 10 -j- 20 K. These results may be compared with the available experimental data on the 
alkali-graphite intercalation compounds (GIC). There is experimental evidence that the critical Temperature 
T c in alkali-GIC C^M (where M is a given alkali metal) grows as x decreases [41]. The number of electrons per 
graphite atom / is related to x by an empirical formula [42] / = n/x where k is the fractional charge transfer. 
A typical value for n is 0.5-^0.6 [43] and hence the filling / = 1.25 corresponds to the alkali-GIC's C2M. 
Under high-pressure, high metal concentration samples such as C&K., C3K, C4Na, CaNa, C2Na, C2Li have 
been synthesized; for C2Na the value of T c is 5 K while for C2Li, T c =1.9 K; in both cases the superconducting 
critical Temperature is of the same order of magnitude of A aS ympt(oo). Quite recently Potassium [35] and 
Lithium [36] have been intercalated also in single- and multi-wall carbon nanotubes up to high concentration 
(the highest metal concentration was obtained with Lithium in C2Li). Our mechanism predicts that the 
binding energy of the W — pairs is bigger in nanotubes than in graphite sheets and this suggests a higher 
critical Temperature for the former, see Fig(5.6). This is also supported by the measurements of a T c w 15 K 
in the 4 Angstrom SWNT by Tang et al. [5]. 



VII. CONCLUSIONS 



Currently, carbon nanotubes superconduct at much lower temperatures than high-T c Cuprates and the two 
kinds of materials are apparently quite different. However, symmetry arguments based on the W = theorem 
tell us that, despite the obvious differences, part of the story must be the same, i.e. by a suitable choice of 
Dirac's characters the on-site Coulomb interaction is utterly turned off. This produces the singlet pairing and 
constrains the ground state spin-orbital symmetry of the interacting system. We presented analytic expressions 
for the effective interaction and obtained the binding energy for (N, N) armchair nanotubes; in the case N = 1 
we verified these predictions numerically and got high-precision agreement. 

Although the results presented in this paper cannot be quantitatively compared with the experimetal data, 
the order of magnitude of the pair binding energy agrees with experiment. Furthermore, the decreasing of the 
binding energy with N is suggested by recent measurements on nanotubes with diameter of few Angstrom [5]. 

The paired state we have obtained here is essentially two-dimensional, that is the transverse direction is 
crucial to have a non-Abelian symmetry group and hence W = pairs; the pairing mechanism uses degenerate 
electronic states that exist in 2d away from half filling. This opens up the interesting possibility that two 
distinct superconducting order parameters appear in the phase diagram, if it turns out that close to half-filling 
there is another one due to a breakdown of the Luttinger liquid. 



14 



REFERENCES 



[1] S. Iijima, Nature (London) 354, 56 (1991). 

[2] M. Bockrath, D. H. Cobden, J. Lu, A. G. Rinzler, R. E. Smalley, L. Balents and P. L. McEuen, Nature 
(London) 397, 598 (1999); Z. Yao, H. W. C. Postma, L. Balents and C. Dekker, Nature (London) 402, 
271 (1999); H. W. C. Postma, M. de Jonge, Z. Yao and C. Dekker, Phys. Rev. B 62, R10653 (2000). 

[3] A. Y. Kasumov, R. Deblock, M. Kodak, B. Reulet, H. Boukiat, I. I. Khodos, Yu. B. Gorbatov, V. T. 
Volkov, C. Journet and M. Burghard, Science 284, 1508 (1999); see also A. Yu. Kasumov, in Proceedigs of 
the School and Workshop on Nanotubes and Nano structures 2000, Santa Margherita di Pula (Ca), Italy. 

[4] A. F. Morpurgo, J. Kong, C. M. Marcus and H. Dai Science 286, 263 (1999). 

[5] Z. K. Fang, Lingyum Zhang, N. Wang, X. X. Zhang, G. H. Wen, G. D. Lee, J. N. Wang, C. F. Chan and 

Ping Sheng, Science 292, 2462 (2001). 
[6] M. Kociak, A. Yu. Kasumov, S. Gue'ron, B. Reulet, I. I. Khodos, Yu. B. Gorbatov, V. F. Volkov, L. 

Vaccarini and H. Bouchiat, Phys. Rev. Lett. 86, 2416 (2001). 
[7] Jeroen W. G. Wildoer, L. C. Venema, A. G. Rinzler, R. E. Smalley and C. Dekker, Nature (London) 391, 

59 (1998); F. W. Odom, J. Huang, P. Kim and C. M. Lieber, Nature (London) 391, 62 (1998). 
[8] N. Hamada, S. Sawada and A. Oschiyama, Phys. Rev. Lett. 68, 1579 (1992). 

[9] X. Blase, Lorin X. Benedict, Eric L. Shirley, and Steven G. Louie, Phys. Rev. Lett. 72, 1878 (1994); C 
L. Kane and E. J. Mele, ibid, 78, 1932 (1997). 
[10] S. J. Fans, M. H. Devoret, H. J. Dai, A. Fhess, R. E. Smalley, L. J. Geerligs and C. Dekker, Nature 
(London) 386, 474 (1997). 

[11] Marc Bockrath, David H. Cobden, Paul L. McEuen, Nasreen G. Chopra, A. Zettl, Andreas Fhess, and 

R. E. Smalley, Science 275, 1922 (1997). 
[12] S. Bellucci and J. Gonzalez, Phys. Rev. B 64, 201106 (2001). 
[13] L. Balents and M. P. A. Fisher, Phys. Rev. B 55, R 11973 (1997). 

[14] Jeno Solyom, in Proceedigs of the School and Workshop on Nanotubes and Nano structures 2000, Santa 

Margherita di Pula (Ca), Italy. 
[15] Y. A. Krotov, D. H. Lee and S. G. Louie, Phys. Rev. Lett. 78, 4245 (1997). 
[16] R. Egger and A. O. Gogolin, Phys. Rev. Lett. 79, 5082 (1997). 
[17] H. Yoshioka and A. A. Odintsov Phys. Rev. Lett. 82, 374 (1999). 
[18] A. Sedeky, L. G. Caron and C. Bourbonnais, Phys. Rev. B 65, 140515 (2002). 
[19] J. Gonzalez, cond-mat/0204171. 

[20] Michele Cini, Gianluca Stelanucci and Adalberto Balzarotti, Eur. Phys. J. BIO, 293 (1999). 

[21] Michele Cini, Enrico Perfetto and Gianluca Stefanucci, Eur. Phys. J. B 20, 91 (2001). 

[22] W. Kohn and J. M. Luttinger, Phys. Rev. Lett. 15, 524 (1965). 

[23] Michele Cini and Adalberto Balzarotti, Phys. Rev. B 56, 14711 (1997). 

[24] Michele Cini and Gianluca Stelanucci, J. Phys.: Condens. Matter 13, 1279 (2001). 

[25] Michele Cini and Gianluca Stelanucci, Solid State Comraun. 117, 451 (2001). 

[26] Gianluca Stelanucci and Michele Cini, J. Phys.: Condens. Matter 14, 2653 (2002). 

[27] Gianluca Stelanucci and Michele Cini, accepted by the Phys. Rev. B. 

[28] Michele Cini, Adalberto Balzarotti, RaHaella Brunetti, Maria Gimelli and Gianluca Stelanucci, Int. J. 

Mod. Phys. B 14, 2994 (2000). 
[29] Michele Cini, Gianluca Stelanucci and A. Balzarotti, Solid State Commun. 109, 229 (1999). 
[30] R. Citro and M. Marinaro, Z. Phys. B 103 153 (1997). 
[31] W. Fettes and I. Morgenstern, Eur. Phys. J. B 9, 9509 (1996). 

[32] J. E. Hirsch, S. Fang, E. Loh and D. J. Scalapino, Phys. Rev. Lett. 60, 1668 (1988). 
[33] C. A. Balseiro, A. G. Rojo, E. R. Gagliano and A. Alascio, Phys. Rev. B 38, 9315 (1988). 
[34] Michele Cini, Adalberto Balzarotti and Gianluca Stelanucci, Eur. Phys. J. B14, 269 (2000). 
[35] M. Bockrath, J. Hone, H. Zettl, P. L. McEuen, A. G. Rinzler and R.E. Smalley, Phys. Rev. B 61, R10606 
(2000). 

[36] B. Gao, A. Kleinhammes, X. P. Fang, C. Bower, L. Fleming, Y. Wu and O. Zhou, Chem. Phys. Lett. 

307, 153 (1999); J. Zhao, A. Buldum, J. Han and J. Ping Lu, Phys. Rev. Lett. 85, 1706 (2000). 
[37] M. P. Lopez Sancho, M. C. Munoz and L. Chico, Phys. Rev. B 63, 165419 (2001). 



15 



[38] A. L. Tchougreef and R. Hoffman, J. Phys. Chem. 96, 8993 (1992). 

[39] M. S. Ferreira, T. G. Dargam, R. B. Muniz and A. Latge', Phys. Rev. B 63, 245111 (2001). 
[40] Lorin X. Benedict, Vincent H. Crespi, Steven G. Louie and Marvin L. Cohen Phys. Rev. B 52, 14935 
(1995). 

[41] J. T. Belash, O. V. Zharikov and A. V. Pal'nichenko, Solid State Comm. 63, 153 (1987); J. T. Belash, A. 
D. Bronnikov, O. V. Zharikov and A. V. Pal'nichenko, Solid State Comm. 64, 1445 (1987); J. T. Belash, 
A. D. Bronnikov, O. V. Zharikov and A. V. Pal'nichenko, Solid State Comm. 69, 921 (1989). 

[42] U. Mizutani, T. Kondow and T. B. Massalski, Phys. Rev. B 17, 3165 (1978). 

[43] R. A. Jishi, M. S. Dresselhaus, Phys. Rev. B 21, 12465 (1992). 

[44] Hannay et al. Phys. Rev. 14, 225 (1965) 



16 



